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1. Introduction 

Recent developments have given us a much better understanding of the degeneracy 
counting formula for ^-BPS dyons in A/" = 4 string compactifications. This formula 
has been considerably refined from its original form in Ref. where it was first 
proposed. One such refinement consists of specifying the integration contour in the 
degeneracy formula and noting that different contours can lead to different answers 
for the degeneracy |^ (for a review, see Ref. [|]). The effect of varying the inte- 
gration contours is in the form of discontinuous jumps in the degeneracy whenever 
the contour crosses a pole in the integration variable and picks up the corresponding 
residue. This has been interpreted as due to the decay of some |-BPS dyons into a 
pair of |-BPS dyons at curves of marginal stability, which are computed using the 
BPS mass formula. 

Because for large charges the decaying states are black holes, a mechanism is 
needed to explain exactly how these decay on curves of marginal stability. The 
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answer turns out to be ^ that |-BPS black holes (for a given set of charges) exist 
both in single-centre and multi-centre varieties. For the latter, the separations of 
the centres are determined by the moduli [0. If we specialise to two-centred dyons 
with both centres being ^-BPS, then it was shown in Ref. that as we approach a 
curve of marginal stability the two centres fly apart to infinity. On the other side 
of the curve the constraint equation has no solutions. This explains (in principle, 
though no method is known to explicitly count states of a two-centred black hole) the 
phenomenon of marginal stability and jumping in the counting formula, in terms of 
the disintegration of two-centred black holes. It should be noted that the degeneracy 
of single-centred black holes with the same charges does not vary across moduli space, 
therefore they exist either everywhere or nowhere. 

In these developments, the only type of marginal decay that plays a role is into 
two |-BPS final states. Also, the only multi-centred black holes needed to complete 
the explanation are those with a pair of |-BPS centres. Though the correspondence 
between these two situations was derived for some special cases, it is believed to 
hold in general, namely for any charge vectors and any point in the entire ^^^y- x 



50(6,22) 



moduli space of A/" = 4 compactifications. 



50(6)x 50(22) 

However, there are many more types of marginal decays in the theory, and in 
one sense they are far more generic. These decays are into a pair of |-BPS final 
states, or into three or more final states each of which can be |-BPS or ^-BPS. In 
another sense these decays are "rare", as it has been shown ||, ||, ^ (at least for 
unit-torsion initial dyons) that they take place on curves of marginal stability that 
have a co-dimension > 1 in the moduli space. Therefore these have been labelled 
"rare decays". In particular they cannot lead to jumps in the degeneracy formula^. 
Nevertheless the existence of such decay modes is of importance in understanding 
the behaviour of dyons as we move around in moduli space, and we will study them 
here for their own sake as well as for possible interesting physical consequences that 
they may turn out to have. 

In Ref. these curves were precisely characterised as circles on the torus moduli 
space. But because of their higher codimension they also need some conditions to 
be satisfied in the remaining part of the moduli space. Though the need for such 
conditions was demonstrated in Refs.|^, ^, the precise conditions have not yet 



^For higher-torsion initial dyons the curves can be of codimension 1, but the degeneracy (or 
rather, index) is still not expected to jump, because of fcrmion zero modes. We will focus largely 
on unit-torsion dyons in this paper. 
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been worked out. In this paper we will completely characterise the codimension > 1 
subspace on which rare decays can take place. 

It is also known that there exist multi-centred dyonic black holes with two |-BPS 
centres, or three or more centres each of which can be |- or |-BPS. However, because 
the degeneracy formula does not jump at curves of marginal stability, these multi- 
centred dyons have not played a role in studies of dyons in jV = 4 compactifications. 
In particular they have not been related to marginal decays into two |-BPS final 
states or multiple final states, and in fact such a relation does not seem necessary 
for the state-counting problem. Nevertheless, in what follows we will argue that the 
relation between curves of marginal stability and multi-centred black holes flying 
apart is quite generic. 

In what follows, we start by briefly reviewing what is known about "rare" 
marginal decays in A/" = 4 compactiflcations. Then we flnd a precise form for the 
constraints on moduli space in order for such rare decays to take place. We exam- 
ine and solve these constraints in a variety of special cases, to give a flavour of what 
they look like. Then using some known results on T-duality orbits, we will obtain the 
constraints in the general case. Next we recursively identify the loci of marginal sta- 
bility for multi-particle decays. Finally we examine the special-geometry formula for 
generic multi-centred black holes and write it in a form that relates their separations 
to curves of marginal stability for n > 2-body decays. 



2. Marginal stability for J\f = 4 dyons 



The electric and magnetic charge vectors of a dyon in an jV = 4 string compactiflca- 
tion are elements of a 28 dimensional integral charge lattice of signature (6, 22). The 
formulae for BPS mass involve a 28 x 28 matrix L, which in our basis will be taken 
to be: 

/ Is ^ 

Ee 
-Ei6, 



(2.1) 



as well as a 28 X 28 matrix M of moduh satisfying MLM^ — L. The inner product 
of charge vectors appearing in the BPS mass is taken with the matrix L + M. In the 
heterotic basis where the compactiflcation is specifled by a constant metric Gij, an 
antisymmetric tensor fleld and constant gauge potentials (where i — 1,2, ■ ■ • ,6 
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and / = 1,2, 



16), this matrix is [11, [T^: 



G' 



l + G-\B + C) 



\ 



\ 



L+M = i + + (G- B + C)G-\G + B + C) (G- B + G)G-^A 

^TQ-i A^G-\G + B + G) A^G-^A 

(2.2) 

Here C is a symmetric 6x6 matrix constructed from A as C = ^A'^A, more concretely 
In this basis we parametrise the charge vectors exphcitly as: 

fo'A fp'\ 



Q 



p 



(6) 

3" 
pill 



pii 

M6) 



(2.3) 



V^(16)/ 



where we have broken up the original vectors into three parts with 6,6 and 16 com- 
ponents respectively. In subsequent discussions we will not explicitly write out the 
subscripts (6), (16) that appear in the above formula. 

The BPS mass formula for |-BPS dyons in A/" = 4 compactifications is as 
follows PI, 111: 



mBPsiQ, Pf 



iQ - fP) o (Q - rP) + 2^2^ ^{Q, P) 



(2.4) 



where 



A(g,F) = (QoQ) (FoF)-(FoQ)2 (2.5) 
The inner products of charge vectors appearing in this formula are of the form: 

Qo P = Q^{L + M)P (2.6) 

The matrix L + M has 22 zero eigenvalues and therefore the inner product only 
contains a projected set of 6 components from the original 28 components of the 
charge vector. Explicitly, the zero eigenvectors take the form: 

(g + b + c A^ \ 

-1 (2.7) 

-1 j 

where each column of the above matrix describes an independent zero eigenvector. 

It is convenient to replace the inner product on charge vectors in Eq. ( p.6| ) by an 
ordinary product acting on some projected vectors. To do this, define L + M as a 



\ 
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^ 

28 X 28 matrix satisfying yL + M yL + M = L + M. This will be ambiguous upto 
a "gauge" freedom but we will select a specific solution that is particularly useful, 
namely: 

/ E-^ E~\G + B + C) E-^A \ 
VL + M = (2.8) 

^ J 

With this matrix it is evident that the projected charges only have their first 6 
components nonzero, namely for any arbitrary vectors Q, P the projected vectors 
Qr, Pr defined by: 

Qr = ^/lTMQ, Pr = VlTMP (2.9) 

are 6-component vectors. The components of these vectors are moduli dependent and 
not quantised. On the projected vectors, one only needs to consider ordinary inner 
products, for example QJ^Qr is equal (by construction) to Q^{L + M)Q. Hence in 
what follows we will denote this quantity either hj QoQ or equivalently by Qr ■ Qr, 
and analogously for other inner products. 

Within the 6-dimensional projected charge space, the electric and magnetic 
charge vectors of the initial dyon span a 2-dimensional plane. Decay of a |-BPS dyon 
into a set of decay products with quantised charge vectors {Q^^\ -P*-^^), ■ ■ ■ , {Q^"'\ P^")) 
can take place only when the plane spanned by the projected charge vectors of each 
decay product coincides with this plane (this is the condition for all states to be 
mutually i-BPS): 

^Qr\ (mi ri\ (Qr 

Si UiJ \Pr^ 

When there are just two decay products and both are |-BPS, the pair of decay 
products defines a 2-plane. Charge conservation then implies that this plane coincides 
with the plane of the original charge vectors, so in this very special case the above 
requirement imposes no conditions on the moduli. Indeed, the numbers mi,ri, Si,ni 
are then integers and the above relation holds between the full (quantised) charge 
vectors, not only the projected ones. Marginal decay takes place on a wall of marginal 
stability whose equation is explicitly known (see Ref.^] and references therein). In all 
other cases, the numbers mi,ri, Si,ni are non-integral and moduli-dependent. Then 
the above condition puts additional constraints on the background moduli M. Our 
goal here is to identify these constraints explicitly. 



(2.10) 



- 5 - 



For a two-body decay into ^-BPS constituents, once the constraints are satisfied 
and we find the numbers mi, ri,si, ni (the corresponding numbers m2, T2, S2, n2 are 
determined by charge conservation) the condition for marginal decay is expressed in 
terms of the curve 0: 

Here we have restricted to the case of unit-torsion dyons, so we have put m = n = 1 
with respect to the notation in Ref. 0. Also, E is defined by: 

E=^(Q«oF_p{i)og) (2.12) 

Interestingly the numerator of this quantity is the Saha angular momentum between 
one of the final-state dyons and the initial state, evaluated with respect to the moduli 
at infinity. Exchanging the role of the two final-state dyons sends E ^ —E. It also 
sends mi — ni (1 — mi) — (1 — rii) = —(mi — ni) and ri, Si — > — ri, —Si. The curve 
of marginal stability is invariant under this set of transformations, as it should be. 

We now turn to the detailed study general two-body decays into |-BPS con- 
stituents. We will find explicit expressions for the numbers mi, ri, si, rti in terms of 
the quantised charge vectors Q, P, Qi, Pi and the moduli M. We will also explicitly 
characterise the loci in moduli space where such rare decays are allowed. 



3. Rare dyon decays 

3.1 Analysis and implicit solution 

It will be useful to define a quartic scalar invariant of four different vectors by: 

A(AS;(7,D)^det( 't°^]=iAoC){BoD)-{AoD)iBoC) (3.1) 

\ B o C B o D I 

As explained above, the "o" product is the mo dull- dependent inner product involving 
the matrix L + M. The above quantity is antisymmetric under exchange of the first 
pair or last pair of vectors, and symmetric under exchange of the two pairs. The 
quartic invariant of two variables defined earlier is a special case of this new invariant: 

A(Q,P) = A(Q,P;Q,P) (3.2) 
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Now start with the following vector equation that is part of Eq. ( 2.10 ): 



Q^R^ = miQR + npR (3.3) 

Contracting this successively with Qr and Pr we find: 

■ Qr = miQl + nQR ■ Pr 
Q^^^ ■ Pr = miQR ■ Pr + nPl (3.4) 

Multiplying the first equation by and the second by Qr ■ Pr and subtracting, we 
find: 

m,A{QR,PR) = A{QR,PR-Q^i,\PR) (3.5) 

which enables us to solve for mi. Repeating this process we can solve for ri,si,ni 
leading to the result: 



V n\ ^ 1 f A{Qr,Pr-Q'^^\Pr) A{Qr,Pr;Qr,Q'^!,'^] 
si nj AiQR,PR)\A{QR,PR;Pl^\PR) A(Qh, P^; Q^, P«) 

It follows that Eq. ( p.lO| ) can be expressed as: 



(3.6) 




1 / A{Qr,Pr;Q^^\Pr) A{Qr,Pr;Qr,Q^'^ 
A{Qr,Pr) \A{QR,PR-p'i\PR) A{Qr,Pr;Qr,P^ 




For fixed, quantised charge vectors Q, P of the initial dyon and Q^^^ , P*-^-* of the first 
decay product (the charge of the second product is determined by charge conserva- 
tion), the above equation provides a set of constraints on the moduli that must be 
satisfied for the \ ^ \ + \ decay to be possible. These constraints together with 
the curve of marginal stability Eq. ( p.ll| ) provide a necessary and sufficient set of 
kinematic conditions for marginal decay. 

In the above form, it is rather difficult to disentangle the constraints or to phys- 
ically understand their significance. Therefore we will consider a number of special 
cases. Along the way we will see the advantage of using T-duality to bring the 
charges into a convenient form and performing the analysis in that basis. Finally we 
write down the explicit constraint equation in the general case, again in the chosen 
T-duality basis. 
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3.2 Explicit solution: special cases 
(i) |-BPS final states 

The case where the decay products are |-BPS should provide no constraints on 
the moduli as this is a "non-rare" decay. This provides a check on our equations. 
Inserting the |-BPS conditions: 



we find that: 



P^'^=hQ^'\ P^^^ = k2Q^'^ (3.8) 



- -^i (3.9, 

Si nij A{Qr,Pr) \k1k2 -ki 



We also have: 





^{Qr, Pr) = ih - k,fA{Q^^\Qf) (3.10) 
Substituting in the above equation, we find: 

f 1 \ 

(3.11) 

At this stage all moduli-dependence has disappeared from the matrix, and equation 
Eq. (|2.10|) indeed provides no constraints on the moduli. Rather, it reduces to an 
identity. It is also easy to see that ki — k2 divides the torsion of the original dyon, 
so if we are also considering the unit-torsion case then ki — k2 = \ and mj, ri, si, rii 
are all manifestly integral [Q. 

(ii) Special charges and moduli 

The next special case we will study has a restricted set of charges. Additionally, 
some of the background moduli are set to a specific value, namely zero in the chosen 
coordinates. We then examine the constraints on the remaining moduli. In choosing 
special values for the moduli, we should in principle avoid loci of enhanced gauge 
symmetry where the dyons we are studying would become massless. 
Let us restrict ourselves to special initial-state charges given by: 

Q' = (g;,o,---,o), q" = (q;',o,---,o), g'" = o (3.12) 

and 

P' = (0, P^, 0, ■ ■ ■ , 0), P" = (0, P^', 0, ■ ■ ■ , 0), P'" = (3.13) 
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Next we set Bij = = as well as Gij = 0, i ^ j- The above restrictions allow 
us to choose the orthonormal frames Eai to be diagonal: 



Eii = Ri, i = 1,2, - ■ ■ ,6 



(3.14) 



with Ri the radii of the six compactified directions in the heterotic basis. 

In the restricted subspace of moduli space that we are considering here, the 



matrix yL + M reduces to: 



VL + M 



v 





E 
























(3.15) 



with E given as in Eq. ( 3.14| ). Therefore the projected initial-state charge vectors 
are: 



+ Q'iR^\ 



Qr 



R 



( \ 

ft + P2R2 





V ) 



\ 



(3.16) 



For this configuration we clearly have Qr ■ Pr = and therefore the quartic invariant 
A is: 

2 / n/ \ 2 

(3.17) 



Ri 



A(g«, Pr) = ( ^ + Q'iRi ) ( £ + ^2 ^2 



We take the decay products to have generic charges Q^^\P^^^ and Q^'^\P^'^^ 
subject of course to the requirement that they add up to Q, P. We then have: 



(1) 



Q'r 



+ Q2 R2 



R2 
.(1)' 



p(i) 

' ^R 



-R^ + n ^1 



(3.18) 
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Now we can compute the quartic invariants appearing in Eq. ( |3.7|) : 



Had we not taken E to be diagonal, the expressions above would have quickly become 
very complicated to write down. 

Inserting the above expressions, and cancelling some common factors, the con- 
straint equation Eq. ( |3.7| ) becomes: 



(3.19) 

These are 6+6 equations. However, the first two components of each set are iden- 
tically satisfied, as one can easily check. This is expected, and follows from the 
structure of Eq. (|2.1CI| ) from which mi,ri,Si,ni were determined. The remaining 
four components of each equation give the desired constraints on the moduli. Be- 
cause of the way we have chosen Q, P, the RHS already vanishes on components 3 
to 6, so the constraint is simply that the LHS vanishes. That in turn sets to zero the 
components 3 to 6 of the vectors and P^^^. Thus we find the constraints: 



(1)' 

+ Ri = 0, z = 3,4,5,6 



Qi , o(i) 



Ri 
p(i)' 

i 

Ri 



+ = 0, 2 = 3,4,5,6 (3.20) 
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If the components of Q'^^\ P^^^ are all nonvanishing, this implies that: 

j = 3,4,5,6 (3.21) 



Q 



(1)' 



p(i)' 



In this special case the constraint equations have some particular features. First of 

all, for generic charge vectors Q^^^ and P^^\ there are no solutions. To have any 

solutions at all, one must choose the charges of the decay products in such a way 

that the second equality in the above equation can be satisfied. In other words, the 
(lY (1)" 

sign of Ql and Ql must be opposite (for i = 3,4,5,6), and the same has to 
be true for P^^\ In this case we find four constraints on the moduli, which fix the 
compactification radii R3, R4^, R^, Rq. 

For this special case, the numbers mi, ri, si, rii in Eq. ( p.lO|) are given by: 



"^1 = —5^ , n = — -, 



Si = rii = (3.22) 

We see that mi, Si depend only on Ri and ri, ni depend only on R2. 

So far the decay products were taken to have generic charges (consistent of course 
with charge conservation). The situation changes if we choose less generic decay 

products. Earlier we took all components of Q^^^ , P^^-* are nonvanishing. However 

(1)' (1)" (1)' (I)" 
if Ql = Ql = P^ = P^ =0 for any i G 3,4,5,6 then the corresponding 



constraint Eq. ( p.20|) is trivially satisfied. In this situation we will have a reduced 



number of constraints. As an example if the above situation holds for all directions 

q(1)' p(l)' 

except 2 = 3 and if — ^ = — fiy^ then there is only a single constraint coming from 



the above equations. The curve of marginal stability provides one more constraint, so 
the decay will take place on a codimension-2 subspace of the restricted moduli space 
in which we are working for this class of examples. The fact that in some situations 
there are no solutions (for example if we do not satisfy that Qi and Q\ have 
opposite signs for i = 3,4,5,6) simply means that our restricted moduli space fails 

to intersect the marginal stability locus in that case. 

(1)' (1)" (1)' (1)" 
If the charges Q] ,Qi , Pi , Pi vanish for alH G 3, 4, 5, 6 then there are no 

constraints (beyond the curve of marginal stability). It is easily seen that this is the 

case where the final states are both |-BPS. 
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(iii) General charges, "diagonal" moduli 



In this subsection we study rare decays allowing for completely general charges Q, P, 
but we will restrict the moduli so that the formulae are tractable. The situation 
turns out to be rather similar to the case studied in the previous subsection. 

Considerable simplification can be brought about in the formulae by using some 
known results on T-duality orbits from Ref. [p!5[ (as reviewed in Appendix A of ||16|| ). 
For this purpose we first change basis from the L matrix used in Ref.||16||: 



L' 



u 





... 


■■ 







0"! 


... 


•• 


• 







... 


■■ 


■ 







ai 


•• 


• 









-Les 





V" 









-LesJ 



(3.23) 



to the one we have defined in Eq. (|2.1|) . Here ai is a Pauli matrix, which occurs 6 
times in the above, and Le^ is the Cartan matrix of Eg. 

In fact using T-duality we will be able to restrict to charge vectors that have the 
last 16 components vanishing, therefore we can ignore these components and work 
in a space of 12-component vectors. We then use a 12 x 12 matrix X that satisfies 



XLX' 



L' 



(3.24) 



to map the equations in Ref.[|T6| to our basis. 

Now the relevant result of T-duality orbits states that any pair of primitive 
charge vectors Q, P can be brought via T-duality to the form: 



g' = (g;,o,---,o), 
p' = (p;,p^,... ,0), 



g" = (g;',o,---,o), q"' = o 
p" = (p;,p^',...,o), P"' = 



(3.25) 



This is close to our previous special case, but with P[, P" turned on. It is no longer 
a special case but represents the general case in a special basis. 

As in the previous example, we restrict the moduli by requiring = Bij = 
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and Gij — 0, i ^ j. Then one finds the projected charges to be: 

f9: + PiRi\ 



Q 



+ Q'iR^\ 




R 



R 



\ / 

The quartic invariant is then found to be: 



Ri 

i + P2R2 



(3.26) 



^R, Pr) - (f + Qi^i) ' (f + P^R^ 



(3.27) 



which is actually the same as in the previous, simpler case where we chose a special 
subset of charges. Computing mi, ri, Si, rii as in the previous subsection and inserting 
them back, the constraint equation can now be written: 



Ri 



Q'[R] 



f + P2R2 ] Q^^ 



(1)' 



Ri 



+ Q?"Ri ] (£ + P^R2 ] - 



R2 



\Ri 



Qr + 



R2 



R 



R 



(1)' 



R2 



+prR2 



Qr + 



^ + QiRi 



P. 



(1)' 



R2 



+ p!,'^"R2 I Pr 



(3.28) 



These equations are slightly more complicated than the previous case for which 
we had Q o P = 0, but the extra complication is only in the first two components, 
which are again trivially satisfied. For the remaining components we find: 



or 

Ri 

p(l)' 
i 

Ri 



+ gP'i?^ = 0, z = 3,4,5,6 



+ P^^^"Ri = 0, 3,4,5,6 



(3.29) 
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These are exactly the same as the constraints we found in the previous case. The 
analysis is therefore also the same: the constraints cannot be satisfied for generic 
charges because our restricted moduli space need not intersect the marginal stability 
locus. When they can be satisfied there are at most four constraints, though there 
will be less if some of the decay product charges vanish. 

3.3 General charges, "triangular" moduli 

In this subsection we restrict the moduli in the most minimal way consistent with 
finding a simple form of the constraint equation. The restriction will be a kind of 
"triangularity" condition: 

{G + B + C)ii = {G + B + C)i2 = 0, i = 3,4,5,6 (3.30) 

with no separate constraints on G, B. A other than the above. 

As before, we use T-duality to put the initial charges into the form of Eq. ( |3.25| ). 
Thereafter, we are still free to make T-duality transformations involving the last four 
components of Q' and Q" and all 16 components of Q'" . The T-duality group is thus 
restricted to an S'O(4,20;2). These transformations will affect the charges of the 
decay products while leaving the initial dyon unchanged. Using them we bring the 
electric charges of the first decay product to the form: 

Q«" = (Qr,Qr,Qr,---,o), 

= (3.31) 

Finally we use an S'0(3, 19; 2) subgroup of T-duality that preserves all the charge 
vectors that we have so far fixed, to bring the magnetic charges of the first decay 
product to the form: 

pW = (pW_p(i)>(.)' p(i)- ..._0), 

p,u'-_(p(..'>(..'>a)''_p(.r_..._o)_ 

P(^)"' = (3.32) 
The charges of the second decay product are determined by charge conservation. 
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Now we use the form of the projection matrix y/L + M and write out Eq. ( p.lO| ) 
exphcitly, after first multiplying through by Eij: 



Qf ^' + {G + B + C),,Q) 



(1)" 



P/'^' + {G + B + C)ijP, 



(!)'■ 



(3.33) 



miQ\ + mi{G + B + C)ijQ'; 

+ nP[ + n{G ^ B + C),jP'^ 
s,Q', + s^{G + B + G),,Q'; 
+ niP^ + ni{G + B + G)i,P^' 

This is a set of 6 + 6 equations. Recall that Gij = AjAj. 

We immediately see that for our choice of T-duality frame for the initial charges, 
as well as using the "triangularity" condition, the RHS of the above equations van- 
ishes for i = 3,4,5,6. Hence we find the constraint equations still in a relatively 
simple form: 

^' + {G + B + G)i,Qf" =0, z = 3, 4, 5, 6 

0, i = 3,4,5,6 (3.34) 



+ (G + P + G)ijPf 



These are then the 4+4 constraints on rare dyon decays, though still with the tri- 
angularity restriction on moduli and in a specific T-duality frame. They must be 
supplemented by the curve of marginal stability, for which we need to know the 
numbers mi , ri , Si , ni . 

The first two components of each line of equations Eq. ( p.33| ) determine the 
values of mi, ri, si, rii. From the first line of those equations we find: 



qW' + (c + P + G)uQ'i'' = miQ[ + m,{G + B + G)uQ- + nP[ + ri(G' + P + G)uPi 



Q 



(1)' 



(G + P + G)2^Qf = nP^ + n{G + P + C)2.Pr 



(3.35) 



Solving for ri from the second equation above, we get: 

g«' + (G + P + C)2.QS'^" 



ri 



P^ + {G + B + G)2iPI' 
Inserting this in the first equation determines mi: 

mi = (^P^ + {G + B + G)2^P:'y\Q[ + iG + B + G)uQ': 



(3.36) 



-1 



qW' + (G + P + G)i,Qf )") (P^ + (G + P + G)2.P'") 
-(qS'^' + (G + P + G)2.g?^") (p; + (G + P + G)i,i^" 



(3.37) 
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Similarly we solve for Si,ni from the second line of Eq. ( 3.33 ) and find: 

pi'^' + (G + s + c)2iP!:^^" 



ni 



Sl 



P^ + {G + B + C)2^Pl' 



P^ + {G + B + C)2.P:'j (Q[ + {G + B + Cl,Q':j x 
(Pf + {G + B + C)uPr) [P2 + iG + B + C)2^Pl') 

- (pi')' + {G + B + G)2^Pi^^") {P[ + r{G + B + G)uP^ 



(3.38) 



Admittedly these are somewhat complicated expressions for the numbers mi, ri, si, rii 
that one needs to plug in to determine the curve of marginal stability on the torus 
moduli space. It is conceivable that a more opportune choice of variables could sim- 
ply them further. Nevertheless, the constraints Eq. ( |3.34| ) on the remaining moduli 
are rather simple. 

3.4 Explicit solution: the general case 

We now turn to the case where the initial and final charges are completely general 
and the moduli are generic as well. Most of the relevant analysis has already been 
done in previous subsections and it only remains to write down the result. However, 
as we will see, the equations rapidly become messy - despite the use of T-duality 
transformations - once we use completely general moduli. 

Let us again start by writing out Eq. (|2.1CI| ) explicitly, but now without any 
condition on the moduli. After multiplying through by Eij, we find the equations: 



+ {G + B + C),,Q) 



P/')' + {G + B 



ChPf 



miQ[ + m,{G + B + G),,Q'; 

+ nP^ + r{G + B + C)ijP;! 
siQ', + s,{G + B + C),,Q'; 

+ niP^ + ni{G + B + G),,P'' 



(3.39) 



which are actually the same as Eq. ( |3.33| ) that we had before. The difference is that 
the RHS no longer vanishes for any of the components (earlier that was guaranteed 
by the triangularity condition that we had assumed on the moduli). Notice that even 
in the most general case, we have gained something by fixing the initial and final 
state charges using T-duality. The last 16 components of these charges have all been 
set to 0, and the result is that most of the terms involving the gauge field moduli 
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have disappeared. The only appearance of these moduh is through Qj = ^i^j 
which in turn only appears in the combination G + B + C. 

This time our strategy will be to choose any 4 equations from the above set of 12 
to determine the variables mi,ni,ri, si. Then in the remaining 8 equations we insert 
these values for the variables and obtain the desired constraint equations. Picking 
the first 2 components for each charge vector, we find: 

qS')' + {G + B + C)i,Qf ^" = miQ; + mi{G + B + C)uQ'l + nP[ + r{G + B + C)uP- 

(3.40) 



gW' + {G + B + C)2iQ\'^" = nP!, + n{G + B + C)2iPl' 



Solving for ri from the second equation above, we get: 

Q^^' + {G + B + C)2^Qf^" 



n 



P^ + (G + B + C)2iPI' 
and inserting this in the first equation, we find mi: 

mi = (P^ + (G + S + C)2iP:') ' {Qi + {G + B + G)uQ':) X 
gS')' + iG + B + C)HQf )") (P^ + (G + S + C)2iP:') 

- (qS^' + {G + B + C)2i gS'^") (Pi' + r{G + B + G)uP^) 
Similarly we solve for Si,ni from the second equation and find: 

pW' + (C + 5 + C)2iP/'^" 



(3.41) 



ni 



P^+{G + B + C)2iPI' 



(3.42) 



(3.43) 



and 



Sl 



(P^ + (G + P + C)2iP:') {Q'i + {G + B + C)uQ'I) X 
(^pW + {G + B + G)hP/')") (p^ + (G + P + C)2iPl') 

- (Pi'^' + (G + P + C)2iPt^") {Pi + {G + B + G)hP,) 



(3.44) 



We feed in these values of mi,ni,ri,Si into the remaining 8 equations to find the 
most general constraint equations on the moduli: 

gf + (G + p + c),,Qf = mi (q: + (G + p + G),,g;') +n(^p: + {G + B + G),,p;) 
p,^')' + (G + p + G),,pf = Sl (g: + (G + p + G),,g;') + m (i^' + (g + p + G),,p;) 

(3.45) 
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here z = 3, 4, 5, 6, and mi, ni, ri, si are given in the above equations. We see that the 
values of mi, ni, ri, si come out the same as in the previous special case, however the 
constraints are much more complicated and - unlike in all the previous special cases 
- depend explicitly on these numbers. Nevertheless, the above equations embody the 
most general kinematic constraints on moduli space to allow a two-body decay of a 
dyon of charges Q, P into |-BPS final state with charges Q*-^-* and P'^^^ (the charges 
of the second state being, as always, determined by charge conservation). It is quite 
conceivable that a more detailed study of possible T-duality bases will allow us to 
further simplify the most general case, and we leave such an investigation for the 
future. 



4. Multi-particle decays 

So far in this work, as well as in previous work^], we have written down conditions for 
decay of a dyon into two |-BPS final states. One could certainly imagine extending 
these considerations to three or more final states. Indeed, it turns out rather simple 
to do so and we will here discuss an iterative way to obtain the relevant formulae. 

Consider the decay of a dyon of charges {Q, P) into n decay products of charges 
(Q(i), P(i)), (Q(2)^ p(2))^ . . . (qW, p(")). The condition for marginality of such a de- 
cay is the condition for the original dyon to go into two decay products of charges 
{Q^^\ -P'-^^) and ^^=2(^*'*^ P^^^)^ along with the condition for the second decay prod- 
uct to further decay into say (Q^^-*, -P*-^-*) and ^^=3(^5''*^ -P''*^)- The latter condition 
must in turn be iterated. Each of these is a two-body decay (with both final states 
being i-BPS) so we already know the condition for each one to take place. The in- 
tersection of all these loci will give the marginal stability locus for the multiparticle 
decay. 

There is a simpler way to iterate the condition. Instead of looking at the curve 
where the second decay product decays into further sub constituents, as above, we 
can simply consider the collection of all marginal stability loci for the decays: 

(t)-(f)*(tf)^ — - 



For each of these, the curve is precisely Eq. ( 2.11|) with the subscript "1" replaced 



by "z". We write it as: 

C{mi,ri,Si,ni) = (^n - ^2g."^' ) (^^ ^) ~ ^ (("^» + + = ^ 

(4.2) 

where 

Ei = ^ [q^'^ o P - pW o (4.3) 

In addition to this curve we have the constraints on the remaining moduh as in Sec. 3 
above. Those too can be expressed in terms of the single decay product labelled "i" . 
Now to find the condition for a multi-dyon decay, we simply take the intersection of 
all these loci of marginal stability. As the number of final states increases, we will 
generically find loci of marginal stability of increasing codimension. 



5. Multi-centred black holes 

It was argued in Refs.[^, |^ that the curves of marginal stability for decays of the 
form: 

i-BPS ^ i-BPS + i-BPS (5.1) 

are also the curves of disintegration for two-centred ^-BPS black holes whose centres 
are individually |-BPS. The method used in these works, which we will summarise 
and extend below, was to use a constraint equation due to Denef [0 to express the 
separation between the centres of such a black hole in terms of charges and moduli. 
Requiring that the separation be infinite places a condition on charges and moduli 
which turns out to be precisely the curve of marginal stability, Eq. (|2.11| ), specialised 
to this decay. 

Now Denef 's constraint equation is not confined to two-centred black holes alone, 
but applies to any number of centres. It has a different limitation: it is defined in 
the context of A/" = 2, rather than Af = 4 compactifications, and relies on special 
geometry. Nevertheless, for the cases to which it applies, we can certainly use it in 
the jV = 4 context. We will do so and will find the perhaps surprising result that the 
curves of marginal stability for generic decays to n final states, which we discussed in 
Sec. 4 above, are precisely reproduced by the constraint equations for multi-centred 
black holes. This suggests a more generic relationship between multi-particle decays 
and multi-centred black holes than has been previously considered. 
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The constraint equation on multi-centred dyons, (see for example Ref. ||^^) reads 
as follows. Let p^'^'^^q^^ be the charges of the i-ih centre where i = 1,2, ■ ■ ■ , N. 
These charges are expressed in the special-geometry basis^. Let the 3- vector fi be 
the location of the z-th centre. And let the moduli be encoded in the standard 
holomorphic special-geometry variables , Fj. Then the constraint equations are: 

Im (F,^) - g« Im (X^) + ^ ^ _ 'i, ^ = (5.2) 

Here the subscript oo indicates that the corresponding moduli are measured at spatial 
infinity (for brevity of notation we will drop it when there is no risk of ambiguity). 
Note that the numerators inside the summation correspond to the Saha angular 
momentum between each pair of centres. 

These are equations for (^) pairwise distances between the centres. We 
analyse them following the procedure in Ref . |^ for the two-centred case. First of all, 
one of the equations is redundant. Adding all the equations, we find: 

p^Im(F,oo)-g/Im(X4) = (5.3) 

where {p^,qi) are the charges of the entire black hole. This provides one real 
constraint on the extra modulus X^. As the above equation is invariant under 
X^ \X' for real X^ , as well as under X^ —X', we see that the magnitude of 
X° is undetermined by this condition, while the phase is determined (in terms of the 
X^,I = 1,2,3) upto a two-fold ambiguity. Another real constraint is now imposed 
in the form of a "gauge condition" : 

X^Fj - X^Fj = -i (5.4) 

This determines the magnitude of X° but leaves intact the two-fold ambiguity in 
the phase. The remaining — 1 equations then provide constraints on the (^) 
separations. 

For the case A^ = 2 we therefore have a single equation, which completely deter- 
mines the separation between the two centres. This works as follows. The relevant 
part of the theory is described by the holomorphic prepotential: 

X^X^X^ , , 
P = (5-5) 

sign in equation (3.2) of Ref.[^ should be corrected so that it reads ^ = — r. This leads to 
some sign changes in other equations there. 

^As we will see, this differs by an interchange of some components from the standard basis used 
in A/" = 4 compactifications. 
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where the are complex scalar fields related to a subset of the K?) x moduli, 
namely r = ri + iT2 describing the 2-torus complex structure, and 

M = diag(i?-^ (5.6) 

describing a 2-parameter subset of the remaining moduli (including the moduli). 
The precise relationship is: 

X^ ^ X^ R , , 



The gauge condition Eq. ( |5.4|) then tells us that: 

= (5.8) 

As in the previous sections, we will consider a dyon with charges {Q,P), but 
now each taken to be 4-component (the first two components should be thought of 
as two of the six Q' and the second two components constitute two of the six Q". 
The charges correspond to unit torsion, namely: 

g.c.d.(Q,F,-P,Q,) = l (5.9) 

We begin by determining the modulus X^ in terms of the T-duality invariants P o 
P, Qo Q, Po Q, where as before the inner products are defined in terms of the moduli 
at infinity, e.g. PoF = P^{L + M)P. 

As promised, we will use the transcription between the natural electric-magnetic 
basis P, Q for the type IIB superstring and the natural basis p^, qi for special geom- 
etry (see for example Ref. |]5|): 

qi = {Qi.Pi.QA.Q2). p' ={P3,-Q3,P2,Pa) (5.10) 
In addition we have: 



while 



Im (Po) = Im (X°r) , Im (Pi) = R^ Im (X°) , 

Im (P2) = ^ Re (XV), Im (P3) = PPRe (X°r) 
R 



Im (X°) = Im (X°) , Im (X^) = - Im (XV) , 

Im (X^) = RR Re (X°) , Im (X^) = ^ Re (X°) 

R 



(5.11) 



(5.12) 
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Inserting these into Eqs.( ^.3|) , (|5.8|) , one finds 



^0 ^ 1 VAf + z(QoPf-QoQ) 
{2V2RT2) VQ^Mbps 



where Mbps is the BPS mass given by Eq. 

Now let us assume our dyon has n centres of charges {Q^"^^ , P*-*-* 

with niijri, Si,ni integers satisfying: 



l,2,---,n (5.14) 



n n n 



'^mi = ^ni = l, ^ri = ^Si = (5.15) 

i=l 1=1 i=l 1=1 

From Eq. ( p.lO|) we find that the charges of the decay products in the qijP"^ basis 
are given by: 



^ (rmQi + ViPi, SiQi + UiPi, niiQi + TiP^, miQ2 + riP^ 



2 



(5.16) 

p^*^^ = {siQs + niP3, -{rriiQs + riP3),SiQ2 + niP2, siQa + UiPi) 



Now the first term in Eq. (p.2|) can be written: 



p^^^ Im {Fi)-qf Im {X^) = RRe (^-X^ X^r 



Si UiJ \ ^ + RP^-t{^ + RP,) 

(5.17) 



The invariants P o P^Q o Q^Q o P are given by: 



Q°Q= yj + PQ3J + + 
Q°^=[^ + ^Q^) + + (^ + ^^^) (5 + ^^^ 

The column vector in Eq. ( ^.17] ) depends on four combinations of Qi, Pi and therefore 
cannot in general be expressed in terms of T-duality invariants. Therefore we restrict 
to the special case, discussed in particular in Ref. for which Qi = Q3 = 0. In this 
case only three independent combinations appear in the column vector and it is easy 
to show that: 



p^'^^ Im (Fi) - qf Im (X^) = PReX° (^-1 r 

Sl^ 



rrii ri\ I y/Q o Q 

QoP+iVA 



Si Tlj 



2V2T2MBPS 



(5.19) 

OljTli^ T-i^ Si^ Tl{j 
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where C{mi,ri, Si,ni) is the curve of marginal stabihty for multiparticle decays, de- 
fined in Eq. gj). 

The numerator of the second term in Eq. (|5.2|), denoted: 



Ji, = p^^'q\'^ -p^^^'qf (5.20) 

is the angular momentum between each pair of decay products evaluated in the 
moduli-independent norm. We will denote the pairwise separation between the cen- 
tres by: 

Lij = \ri - fj\ (5.21) 

Note that J'ij = —J^i and Lj^ = L^j. 

Inserting the above results into Eq. ( p. 21 ), one finds that it can be expressed as 
follows: 

C^ + Y.f^ = ^ (5.22) 

where 



Ci = J^' ^ — C(mi, n, Si, Ui) (5.23) 
V2 T2 Mbps 

Clearly the first term in Eq. (|5.22|) depends only on the charges of a single centre 
(as well as the initial charges) while the second term depends on the charges of 
a pair of centres. Note that we have X]?,^* ~ 0- Thus we have shown that the 
curves of marginal stability for multi-centred decays appear also from considerations 
of multi-centred black holes and the constraints on the locations of their centres. 

In the special case considered previously [§, ^ where the dyon has two |-BPS 
centres, the corresponding curve of marginal stability is of codimension 1. In this 
case it is known that the degeneracy of states jumps as we cross the curve. From 
the supergravity point of view, it was suggested in the Af = 2 context in Ref . ^ and 
shown more explicitly in the present A/" = 4 context in Refs.|^, that this decay 
occurs as a result of the two centres flying apart to infinity at a curve of marginal 
stability. This is seen by specialising Eq. ( p. 221 ) to this case. As long as J7i2 7^ 0, the 



separation L12 — * 00 when Cj — > 0. Moreover for a fixed sign of J'12, the separation 
L12 can be positive only on one side of the curve of marginal stability. On the other 
side it is negative, which indicates that the corresponding two-centred black hole 
does not exist. 

Now let us return to the more general case where there are two centres but both 
are |-BPS. As we have seen, in this case the locus of marginal stability is not a wall 
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in moduli space, but rather a curve of codimension > 2. Therefore the degeneracy 
formula cannot jump as one crosses the curve. Hence one need not have expected 
any relationship between marginal decays and multi-centred dyons. Nevertheless, we 
see that Eq. ( |5.22| ) continues to hold in the more general case (with the limitation 
that the charges are those that can be embedded in an A/" = 2 compactification). 

We interpret this as evidence that the relationship between dyon decay and the 
disintegration of multi-centred black holes holds more generally than required by 
the degeneracy formula. Therefore we conjecture that even with the most general 
charges, n-centred black holes exist in A/" = 4 string compactifications with generic 
|-BPS centres for which Eq. ( ^.22|) holds true. It would be worth trying to prove that 
this is the case, or else to show that such solutions do not exist beyond the cases that 
can be embedded in the charge space and moduli space of J\f = 2. An intermediate 
possibility also exists: that in A/" = 4 compactifications such multi-centred black 
holes do exist with arbitrary charges, but only on a subspace of the moduli space. 

Examining Eq. ( p.22|) one sees that if the marginal stability condition = is 
satisfied for a particular i, then we must have: 

E I" = (5-24) 

One possible solution is to have Lij oo for all j ^ i. This means the ith centre 
has been taken infinitely far away from all the others, in agreement with the picture 
of marginal decay that we developed in Section 4 above. Since the pairwise Saha 
angular momenta J'ij = P*-*^ • Q^^^ — P*--'-' ■ Q*-*^ cannot all be positive in every equation 
(since J'ij = —Jji) there could be other configurations where the C, = 0, except in 
the case of two centres. It is not clear to us how these other solutions should be 
interpreted. 

Note that in the above equation the angular momentum is measured with respect 
to moduli-independent inner product P ■ Q = P^LQ unlike the angular momentum 
appearing in the curve of marginal stability Eq. (|4.2|) which is computed using the 
mo dull- dependent inner product PoQ = P'^[L + M)Q. One may think of the latter 
evolving to the former as we follow the attractor flow from infinity to the horizon of 
the black hole. However it would be nice to have a better physical understanding of 
the role of dyonic angular momenta in these discussions^. 

^As is well-known, the dyonic angular momentum plays a physical role in the wall-crossing 
formulae |l^, |, |, |l§ that describe how the de generacy jumps, but in the present discussions 
there are no walls or jumps. 
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6. Conclusions 

In this work we have obtained the loci of marginal stability for decays of ^-BPS dyons 
into any number of BPS constituents in AT = 4 string compactifications. These loci 
appear as equations constraining the 132 + 2 moduli, more precisely as a curve of 
marginal stability in the upper-half-plane that represents a torus moduli space (in 
the basis of type IIB on x T^, this is the geometric torus) as well as some more 
complicated equations on the remaining moduli. While in this paper we worked with 
unit-torsion initial dyons, it should be quite straightforward to extend our results to 
general torsion. We showed how to extend our analysis to multi-particle decays, and 
found a relation between the loci of marginal stability obtained in this way and the 
supergravity constraints on pairwise separations of the centres of multi-centred black 
holes. 

The physical role of "rare" marginal dyon decays, namely all those other than of 
a |-BPS dyon into two |-BPS dyons, has yet to be explored. Because such decays 
take place on loci of codimension > 2 in moduli space, they do not form "domain 
walls" across which the degeneracy can jump. Therefore they do not affect the basic 
entropy or dyon counting formulae. However it is certainly possible that they have 
other interesting physical effects which may emerge on further investigation. 
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